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The consensus of opinion in cosmology is that the Universe is currently undergoing a period of 
accelerated expansion. With current and proposed high precision experiments it offers the hope of 
being able to discriminate between the two competing models that are being suggested to explain 
the observations, namely a cosmological constant or a time dependent 'Quintessence' model. The 
latter suffers from a plethora of scalar field potentials all leading to similar late time behaviour 
of the universe, hence to a lack of predictability. In this paper, we develop a model independent 
approach which simply involves parameterizing the dark energy equation of state in terms of known 
observables. This allows us to analyse the impact dark energy has had on cosmology without the 
need to refer to particular scalar field models and opens up the possibility that future experiments 
will be able to constrain the dark energy equation of state in a model independent manner. 

PACS numbers: 98.70.Vc,98.80.Cq 



INTRODUCTION 



Current cosmological observations suggest that the en- 
ergy density of the universe is dominated by an unknown 
type of matter, called 'dark energy', Q 0, H 0, S @. 
If correct, it is characterized by a negative value for 
the equation of state parameter w (close to —1) and 
is responsible for the present accelerating expansion of 
the Universe. The simplest dark energy candidate is 
the cosmological constant, but in alternative scenarios 
the accelerated phase is driven b y a dynamical scalar 
field called 'Quintessence' [li @, Although recent 

analysis of the data provides no evidence for the need 
of a quintessential type contribution 0, 0, 0], the 
nature of the dark energy remains illusive. In fact the 
present value of its equation of state wq is constrained 
to be close to the cosmological constant one, but the 
possibility of a time dependence of w q o r a coupling 
with cold dark matter for example [1J| cannot be 
excluded. Recent studies have analysed our ability 
to estimate wq with high redshift object observations 
[TH Hil IrR Ha . Jis! and to reconstruct both the time 
evolution of wq YM. Ell El IH El El Er^ as well as the 
scalar field potential V{Q) |27l l28l l2flj. An alternative 
method for distinguishing different forms of dark energy 
has been introduced in [3(j , but generally it appears that 
it will be difficult to really detect such time variations 
of wq even with the proposed SNAP satellite |3ll |32| . 
One of the problems that has been pointed out is that, 
previously unphysical fitting functions for the luminosity 
distance have been used, making it difficult to accurately 
reproduce the properties of a given quintessence model 
from a simulated data sample [2fl Erj l34| . A more 
efficient approach consists of using a time expansion 
of wq at low redshifts. For instance in [2ll EflL Eri | a 
polynomial fit in redshift space z was proposed, while 



in [33l 13 1| a logarithmic expansion in z was proposed 
to take into account a class of quintessence models with 
slow variation in the equation of state. However even 
these two expansions are limited, in that they can not 
describe models with a rapid variation in the equation 
of state, and the polynomial expansion introduces a 
number of unphysical parameters whose value is not 
directly related to the properties of a dark energy 
component. The consequence is that their application 
is limited to low redshift measurements and cannot be 
extended for example to the analysis of CMB data. An 
interesting alternative to the fitting expansion approach, 
has recently been proposed [$3, in which the time 
behavior of the equation of state can be reconstructed 
from cosmological distance measurements without 
assuming the form of its parametrization. In spite of 
the efficiency of such an approach, it does not take into 
account the effects of the possible clustering properties of 
dark energy which become manifest at higher redshifts. 
Hence its application has to be limited to the effects 
dark energy can produce on the expansion rate of the 
universe at low redshifts. On the other hand, it has been 
argued that dark energy does not leave a detectable 
imprint at higher redshifts, since it has only recently 
become the dominant component of the universe. Such 
a statement, however, is model dependent, on the face 
of it there is no reason why the dark energy should 
be negligible deep in the matter dominated era. For 
instance CMB observations constrain the dark energy 
density at decoupling to be less then 30 per cent of 
the critical one 36]. Such a non negligible contribution 
can be realized in a large class of models and therefore 
cannot be a priori excluded. Consequently it is of crucial 
importance to find an appropriate parametrization for 
the dark energy equation of state that allows us to 
take into account the full impact dark energy has on 
different types of cosmological observations. In this 



2 



paper we attempt to address this very issue by providing 
a formula to describe wq, that can accommodate 
most of the proposed dark energy models in terms of 
physically motivated parameters. We will then be in a 
position to benefit from the fact that the cosmological 
constraints on these parameters will allow us to infer 
general information about the properties of dark energy. 



DARK ENERGY EQUATION OF STATE 

The dynamics of the quintessence field Q for a given 
potential V is described by the system of equations: 



Q + 3HQ + ^- = 0, 
dQ 



and 



H 2 = 



8vrG 



O 2 

Pm +Pr + -Y + V(Q) 



(1) 



(2) 



where p m and p r are the matter and radiation energy 
densities respectively and the dot is the time derivative. 
The specific evolution of wq(cl), where a is the scale fac- 
tor, depends on the shape of the potential, however there 
are some common features in its behaviour that can be 
described in a model independent manner and which al- 
low us to introduce some physical parameters. As a first 
approach, we notice that a large number of quintessence 
models are characterized by the existence of the so called 
'tracker regime'. It consists of a period during which the 
scalar field, while it is approaching a late time attrac- 
tor, evolves with an almost constant equation of state 
whose value can track that of the background component. 
The necessary conditions for the existence of tracker so- 
lutions arising from scalar field potentials has been stud- 
ied in 0, 1371 . In this paper, we consider a broad class 
of tracking potentials. These include models for which 
wg (a) evolves smoothly, as with the inverse power law 
|38j . V(Q) ~ l/Q" (INV) and the supergravity inspired 
potential [U, V{Q) ~ 1/Q a e^/ 2 (SUGRA). Late time 
rapidly varying equation of states arise in potentials with 
two exponential functions |io|. V ~ e~ aQ + e /3Q (2EXP), 
in the so called Albrecht & Skordis' mode l lil (AS) and 
in the model proposed by Copeland et al. 42] (CNR). To 
show this in more detail, in fig.l we plot the equation of 
state obtained by solving numerically Eq. (JIJ and Eq. (J2J 
for each of these potentials. 

There are some generic features that appear to be 
present, and which we can make use of in our attempts to 
parameterize wq . For a large range of initial conditions 
of the quintessence field, the tracking phase starts be- 
fore matter-radiation equality. In such a scenario wq (a) 
has three distinct phases, separated by two 'phase tran- 
sitions'. Deep into both the radiation and matter dom- 
inated eras the equation of state, wq(cl), takes on the 
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FIG. 1: Evolution of wq against the scale factor for an in- 
verse power law model (solid blue line), SUGRA model (|3^|) 
(dash red line), two exponential potential model ([4Q|) (solid 
magenta line), AS model f|4l|1 (solid green line) and CNR 
model ( 42] ) (dot orange line). 
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FIG. 2: Plot of ii)g(a) best fit for different potentials 



values Wq and Wq respectively, values that are related 
to the equation of state of the background component 
w B through HE3]: 



w B - 2(r - i) 
2r-i 



(3) 



where Y = V"V/(V') 2 and V = dV/dQ etc. For the case 
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of an exponential potential, r = 1, with wq — wb, but 
in general wq ^ wb ■ Therefore if we do not specify the 
quintessence potential the values of Wq and Wq should 
be considered as free parameters. 

The two transition phases can each be described by 
two parameters; the value of the scale factor a r c ,m when 
the equation of state wq begins to change and the width 
A-r,m °f the transition. Since T is constant or slowly vary- 
ing during the tracker regime, the transition from Wq to 
Wq is always smooth and is similar for all the models (see 
fig. 1). To be more precise, we have found that a r c ~ 10~ D 
and A r ~ 10~ 4 during this transition, the former num- 
ber expected from the time of matter-radiation equality 
and the latter from the transition period from radiation 
to matter domination. However, when considering the 
transition in wq from Wq to the present day value Wq, 
we see from fig. 1 that this can be slow (0 < a™/A m < 1) 
or rapid (a™/A m > 1) according to the slope of the 
quintessence potential. For instance in models with a 
steep slope followed by a flat region or by a minimum, 
as in the case of the two exponentials, the AS poten- 
tial or the CNR model, the scalar field evolves towards a 
solution that approaches the late time attractor, finally 
deviating from the tracking regime with the parameter T 
rapidly varying. In contrast the inverse power law poten- 
tial always has a slower transition since T is constant for 



where oq = 1, and the value of and a r and be ar- 

bitrarily chosen in the radiation and matter era because 
of the almost constant nature of wq during those eras. 
For example in our simulations we assumed a r = 1CP 5 
and a, n = 10~ 3 . In table 1 we present the best fit pa- 
rameters obtained minimizing a chi-square for the differ- 
ent models of fig.l, and in fig. 2 we plot the associated 
functions Wq (a) . It is encouraging to see how accurately 
the Fermi-Dirac functions mimic the exact time behavior 
of wq(ci) for the majority of the potentials. In fig. 3 we 
plot the absolute value of the difference Aiu(a) between 
wq(o) and Wq(o). The discrepancy is less then 1% for 
redshifts z < 10 where the energy density of the dark 
energy can produce observable effects in these class of 



all times. Given these general features we conclude that 
the behavior of wq (a) can be explained in terms of func- 
tions, Wg(a), involving the following parameters: Wq, 
Wq, Wq, a™ and A m . The authors of |43j have recently 
used an expansion in terms of a Fermi-Dirac function in 
order to constrain a class of dark energy models with 
rapid late time transitions. In what follows we find that 
a generalisation of this involving a linear combination of 
such functions allows for a wider range of models to be in- 
vestigated. To be more precise, we propose the following 
formula for Wq (a) : 



w p Q (a) = F 1 j r {a) + F 2 f m (a) + F 3 , (4) 

with 

fr,m{a) = a _ tt r,m ■ (5) 

1 + e 

The coefficients F%, F 2 and F3 are determined by de- 
manding that Wq(o) takes on the respective values Wq, 
Wq, Wq during radiation (a r ) and matter (a m ) domina- 
tion as well as today (ao). Solving the algebraic equations 
that follow we have: 



(6) 



(7) 



(8) 
I 

models and it remains below 9% between decoupling and 
matter-radiation equality. Only the CNR case is not ac- 
curately described by Wq{a) due to the high frequency 
oscillations of the scalar field which occur at low rcdshift 
as it fluctuates around the minimum of its potential. In 
fact these oscillations are not detectable, rather it is the 
time-average of wq (a) which is seen in the cosmological 
observables, and can be described by the corresponding 
Wq(a). There are a number of impressive features that 
can be associated with the use of Wg(a) in Eq. J1J. For 
instance it can reproduce not only the behavior of models 
characterized by the tracker regime but also more general 
ones. As an example of this in fig. 4 we plot Wq(o) cor- 



(itffi - W Q ) (f m (ao) - /m(flr)) - {Wq - W r Q ) (f m (a m ) ~ fmM) 

(fr(a m ) - fr{a r )) (fm(ao) - f m (a r )) - (fr(ao) - fr(a r )) (fm{a m ) - f m (a r )) ' 
p 2 = <~ W Q _ Fi Mao) -Mar) | 

fm{ao) — }m{a r ) M( a 0) — fmiflr) 

F 3 = w r Q- FiM a r) - F 2 f m {a r ), 
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responding to three cases: a K-essence model 44] (blue 
solid line); a rapid late time transition jl^ (red dash-dot 
line) and finally one with an equation of state Wq < — 1 
(green dash line). The the observational constraints on 
Wq, Wq, Wq, a™ and A m lead to constraints on a large 
number of dark energy models, but at the same time it 
provides us with model independent information on the 
evolution of the dark energy. It could be argued that the 
five dimensional parameter space we have introduced is 
too large to be reliably constrained. Fortunately this can 
be further reduced without losing any of the essential de- 
tails arising out of tracker solutions in these Quintessence 
models. In fact nucleosynthesis places tight constraints 
on the allowed energy density of any dark energy compo- 
nent, general ly fo rcing them to be negligible in the radi- 
ation era |47l l48j . The real impact of dark energy occurs 
after matter-radiation equality, so we can set Wq — Wq. 
Consequently we end up with four parameters: Wq, Wq, 
a™ and A m . Although they increase the already large 
parameter space of cosmology, they are necessary if we 
are to answer fundamental questions about the nature of 
the dark energy. The parameters make sense, if wq{cl) 
evolves in time, we need to know when it changed (a™), 
how rapidly (A m ) and what its value was when it changed 
(wq ). Neglecting the effects during the radiation dom- 
inated era it proves useful to provide a shorter version 
of Eq. (@J, in fact since we can neglect the transition 
from radiation to matter dominated eras, then the linear 
combination Eq. Q can be rewritten as *: 



TABLE I: Best fit values of the parameters of the expansion 









T 


a? 


A m 


INV 


-0.40 


-0.27 


-0.02 


0.18 


0.5 


SUGRA 


-0.82 


-0.18 


0.10 


0.1 


0.7 


2EXP 


-1 


0.01 


0.31 


0.19 


0.043 


AS 


-0.96 


-0.01 


0.31 


0.53 


0.13 


CNR 


-1.0 


0.1 


0.32 


0.15 


0.016 
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FIG. 3: Absolute value of the difference between wq(cl) and 
TOq(a) for the models of fig.l. 



<?( a ) = w q + K> - w q) x x -; — 5^- 

l + e szr 1-e^ 

(9) 

As we can see in fig. 5, the relative difference between 
the exact solution wq(cl) of the Klein-Gordon equation 
and Eq. © is smaller than 5% for redshifts z < 1000, 
therefore it provides a very good approximation for the 
evolution of the quintessence equation of state. Both 
Eq. (JIJ and Eq. @ are very useful in that they allow 
us to take into account the clustering properties of dark 
energy (see for instance |46j ) and to combine low redshift 
measurements with large scale structure and CMB data. 
We would like to point out that one of the key results in 
this paper is the fact that our expansion for the equation 
of state allows for a broader range of redshift experiments 
than has previously been proposed. 



* We thank Eric Linder for pointing this out to us. 



CONCLUSIONS 

The evidence for a present day accelerating universe 
appears to be mounting, and accompanied with it is the 
need to understand the nature of the dark energy that 
many believe to be responsible for this phenomenon. Of 
the two possibilities proposed to date, a cosmological con- 
stant or a Quintessence scalar field, the latter suffers in 
that there are a plethora of models that have been pro- 
posed, all of which satisfy the late time features, that of 
an accelerating universe. Yet there is no definitive parti- 
cle physics inspired model for the dark energy therefore 
an other route should be explored, trying to determine 
the equation of state that the dark energy satisfies in a 
model independent manner. If this were to work it would 
allow us to discuss the impact dark energy has had on 
cosmology without the need to refer to particular dark 
energy scenario. In this paper, we have begun addressing 
this approach. We have introduced a parameterization of 
the dark energy equation of state Wq(o) which involves 
five parameters and shown how well they reproduce a 
wide range of dark energy models. By estimating their 
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FIG. 4: Time evolution of Wq (a) as in the case of K-essence 
(blue solid line), late time transition (red dash-dot line) end 
with Wq < — 1 (green dash line). 
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FIG. 5: Absolute value of the difference between wq(<i) and 
the low redshift formula Eq. J5J for the models of fig.l. 



value from cosmological data we can constrain the dark 
energy in a model independent way. This could be impor- 
tant in future years when high precision CMB and large 
scale structure observations begin to probe medium to 
large redshifts, regions where differences in the features 
of the Quintessence models begin to emerge. Although 
the low redshift data indicate that it is currently im- 
possible to discriminate a cosmological constant from a 



Quintessence model, moving to higher redshifts may re- 
sult in evidence for a time varying equation of state in 
which case, we need to be in a position to determine such 
an equation of state in a model independent manner. We 
believe this paper helps in this goal. 

We would like to thank Carlo Ungarelli, Nelson J. 
Nunes, Michael Malquarti and Sam Leach for useful dis- 
cussions. In particular we are grateful to Eric Linder for 
the useful suggestions he made. PSC is supported by a 
University of Sussex bursary. 



[1 
[2 
[3 
[4 
[5 
[6 

ir 

[8 
[9 
[10 

[11 

[12 

[13 
[14 
[15 
[16 

[17; 
[is; 

[19 

[20 

[21 
[22 
[23 
[24 

[25 

[26 
[27 
[28 
[29 

[30 
[31 

[32 
[33 
[34 



Perlmutter, S. et al. 1999, Astrophys. J. , 517, 565 
Riess, A. et al. 1999, Astrophys. J. , 117, 707 
De Bernardis, P. et al. 2000, Nature 404, 955-959 
Netterfield, C. B. et al. 2001, jasiro-ph/0104460] 
Pryke, C. et al. 2001, Ast rophys. J. , 568, 46 
Efstathiou, G. et al. 2001, |astro-ph /0109152 
Wetterich, C. 1988, Nucl. Phys. B 302, 645 
Ratra, B. and Pebbles, P.J.E. 1988, Phys. Rev. D 37, 
3406 

Caldwell, R. R., Dave, R., and Steinhardt, P. J. 1998, 
Phys. Rev. Lett. 80, 1582-1585 

Ferreira, P. G, and Joyce, M. 1998, Phys. Rev. D58, 
023503 

Corasaniti, P. S. and Copeland, E. J. 2002, Phys. Rev. 
D65, 043004 

Bean, R. and Melchiorri, A. 2002, Phys. Rev. D65, 
041302 

Hannestad, S. and Morts ell, S. 2002, astro-p h/020509 
Amendola, L. et al. 2002, astro-ph/0205097 
Waga, I. and Miceli, A. P. 1999, Phys. Rev. D59, 103507 
Newman, J. A. and Davis, M. 2000, Astrophys. J. Lett., 
534, Lll 

Newman, J. A. and Davis, M. 2002, Astrophys. J. , 564, 
567 

Lima, J. A. S. and Alcaniz, J. S. 2002, Astrophys. J. , 
566, 15 

Calvao, M. O., de Mello Neto, J. R. T. and Waga, I. 
2002, Phys. Rev. Lett. 88, 9 

Huterer, D. and Turner, M. 2001, Phys. Rev. D64, 123527 

Astier, P. 2001, Phys. Lett. B, 500,8 

Barger, V. and Marfatia, D. 2001, Phys. Lett. B, 498, 67 

Goliath, M. et al. 2001, astro-ph/0104009 

Maor, I., Brustein, R. and Steinhardt, P. 2001, Phys. 

Rev. Lett, 86, 6 

Weller, J. and Albrecht, A. 2001, Phys. Rev. Lett. 86, 
1939 

Weller, J. and Albrecht, A. 2002, Phys. Rev. D65, 103512 
Saini, T.D. et al. 2000, Phys. Rev. Lett. 85, 1162 
Huterer, D. and Turner, M. 1999, Phys. Rev. D60, 081301 
Chiba, T. and Nakamura, T. 2000, Phys. Rev. D62, 

121301 

Sahni, V. et al. 2002, astro-ph/0201498 

Kujat, J. et al. 2001, |astro-ph/0112221*f Astrophys. J. in 

press 

Maor, I. et al. 2001, astro-ph/0112526 

Efstathiou, G. 2000, MNRAS 342, 810* 

Gerke, B. F. and Efstathiou, G. 2002, |astro-p h/0201336 
MNRAS in press 



6 



[35] Huterer, D. and Starkman, G. 2002, astro-ph/0207517 
[36] Bean, R., Hansen, S. H. and Melchiorri, A. 2001, 

Phys.Rev. D64 103508 
[37] Steinhardt, P. J., Wang, L. and Zlatev, I. 1999, Phys. 

Rev. D59, 123504 
[38] Zlatev, I., Wang, L., and Steinhardt, P. J. 1999, Phys. 

Rev. Lett. 82, 896-899 
[39] Brax, P. and Martin, J. 1999, Phys. Lett. B468, 40-45 
[40] Barreiro, T., Copeland, E. J., and Nunes, N. J. 2000, 

Phys. Rev. D61, 127301 
[41] Albrecht, A., and Skordis, C. 1999, Phys. Rev. Lett. 84, 

2076-2079 



[42] Copeland, E. J., Nunes, N. J., and Rosati, F. 2000, Phys. 

Rev. D62, 123503 
[43] Bassett, B. et al. 2002, astro-ph/0203383 
[44] C. Armendariz-Picon, V. Mukhanov and P. J. Steinhardt, 

2000, Phys. Rev. Lett. 85, 4438 
[45] L. Parker and A. Raval, Phys. Rev. D60 123502 (1999) 
[46] R. Dave, R. R. Caldwell, P. J. Steinhardt, Phys.Rev. D66 

(2002) 023516 

[47] Bean, R., Hansen, S. H. and Melchiorri, A. 2001, 

Phys.Rev. D64, 103508 
[48] Yahiro, M. et al. 2002, Phys.Rev. D65, 063502 



